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> 1 INTRODUCTION 



ABSTRACT 

We compare orbits in a thin axisymmetric disc potential in MOND to those in a thin 
disc plus near-spherical dark matter halo predicted by a ACDM cosmology. Remark- 
ably, the amount of orbital precession in MOND is nearly identical to that which 
occurs in a mildly oblate CDM Galactic halo (potential flattening q=0.9), consistent 
with recent constraints from the Sagittarius stream. Since very flattened mass distri- 
butions in MOND produce rounder potentials than in standard Newtonian mechanics, 
we show that it will be very difhcult to use the tidal debris from streams to distinguish 
between a MOND galaxy and a standard CDM galaxy with a mildly oblate halo. 

If a galaxy can be found with either a prolate halo, or one which is more oblate 
than (7 ~ 0.9 this would rule out MOND as a viable theory. Improved data from the 
leading arm of the Sagittarius dwarf - which samples the Galactic potential at large 
radii - could rule out MOND if the orbital pole precession can be determined to an 
accuracy of the order of ±1°. 
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Ever since Zwicky's seminal work in the 1930's it has been 
known that there is a disparity between the mass of galaxies 
as measured dynamically and the mass inferred from the vis- 
ible light. The standard explanation for this missing matter 
is to invoke one or many weakly interacting massive particles 
which form early in the universe and to first order interact 
only via gravity (see e.g. Bcrgstrbm (1998)). This is known 
as cold dark matter or CDM theory. However, since none 
of the well motivated candidate particles has yet been de- 
tected, it is important to also consider alternative theories 
as an explanation for the missing matter. 

One such alternative theory, is a modification of stan- 
dard Newtonian gravity (MOND) for acceleratio ns below 
some char acteristic scale, g p ~ 1.2 x 10~^''ms~^ iMilgromI 
Jl983h and lMcGauehl J^ioi)). MOND was first suggested by 
Milgrom in 1983 as a modified inertia theory, but since then 
has be en expanded into a self cons istent Lagrangian field 
theory jBekenstein fc Milgromll984l) and, more recently, has 
been place d on a firm footin g within the context of general 
relativity jBekeri steinl 12004) . This last point is of particu- 
lar interest since Bekcnstoin (2004^) has managed to address 
many of the conceptual problems that have plagued MOND 
over the past two decades and shown that the theory can 
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be consistent with gravitational lensing and other general- 
relativistic phenomena. 

In this paper we compare the potential of the Milky 
Way as predicted by MOND and CDM models. In the for- 
mer the potential arises from a fiattened disk of baryons 
whereas the latter potential is primarily from an extended 
spheroidal distribution of dark matter. A perfectly spheri- 
cal potential has orbits wh ich are confined to lie on planes 
jBinnev fc Tremainelll987^ . By contrast, orbits in axisym- 
metric potentials generally show precession of their orbital 
planes (this will be true for all orbits which are not exactly 
planar or exactly polar). 

Ilbata et alJ i200ll) and m ore r ecently iMaiew ski et alJ 
J2003I) . I Johnston et all 12005") and iLaw et alj ||200^ have 
studied the tidal debris from the Sagittarius dwarf galaxy 
and calculated hkely orbits for the galaxy and its stellar 
debris. They find that the precession of the orbital plane is 
small (~ 10°) and is consistent with a galactic halo potential 
which is only mildly oblate (g = 0.9—0.95) (although see also 
iHelmi (200_4) and section|51in this paper). In MOND, where 
all of the gravity comes only from the disc, the potential may 
be much flatter leading to far more precession than is ob- 
served. In this way, tidal debris from infalling satellites such 
as the Sagittarius dwarf could provide strong constraints on 
any altered gravity theory which supposes that all gravity 
is produced only by the visible light. 

This paper is organised as follows: In section |5| we dis- 
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cuss the Galactic MOND potential. In section |3 we outline 
the initial conditions and orbit solver. We use four models: 
the CDM model which has a spherical dark matter halo, the 
f095CDM and f09CDM models which have shghtly oblate 
halos {q = 0.95 and q = 0.9) and the MOND model, where 
all of the gravitational potential comes from the disc. In sec- 
tion|l|we present our results for two different orbits: the first 
is motivated by the orbit of the Sagittarius dwarf and the 
second is a small-pericentre orbit chosen to sample a wide 
range of the Galactic potential. In section |K| we briefly dis- 
cuss the significance of these results and relate our work to 
previously published studies. Finally, in section El we present 
our conclusions. 



2 THE MOND POTENTIAL 

The MOND field equations lead to a modified, 
non- linear, version of Po isson's equation given by 
jBekenstein fc MUgromlll984) : 



V • [M(|V($)|/ao)V$] = AtvGp 



(1) 



Where p is the density, $ is the scalar field for MON- 
Dian gravity and ao is the acceleration scale below 
which gravity deviates from standard Newtonian behaviour. 
The unknown function V(cl?)|/ao) parameterises the 
change from Newtonian to MONDian gravity and is usu- 
ally given ph enomcnologically by p,{x) = x{l + x^)~^^'^ 
jBekenstein fc Mikrom 19841 . 

Equation Q is in general extremely difficult to solve, 
not least because it is trivial to show that substituting 
$ — > $1 + $2 does not give p pi + p2- This means 
that solutions cannot be superposed as in normal Newto- 
nian mechanics. Every mass configuration will have its own 
unique potential which should be determined by (numeri- 
cally) inverting equation ^ Thus, while some authors (see 
e.g. ,Kncbc & Gibson (2004)) have made valiant efforts to 
adapt N-body integrators to work in MOND, these can only 
be, at best, an approximation. 

However, eq uationUmay be solved in extremely special 
cases. Following iBrada fc MilgromI 1^9^, notice that we 
may write the MONDian gravitational field {g = V$) as 
the sum of the Newtonian gravitation field (gjv — V"I>]v) 
and a curl field: 



(2) 



The curl field will trivially vanish for planar, spherical 
or cylindrical symmetry giving (in exact a greement with t he 
modified inertia interpretation of MOND lMilgrom|[l983) '): 

M(ip!/ao)ff = ffjv (3) 

)~^''^ as above 

(4) 



Which, substituting for p{x) = x(l 
and inverting gives: 

,1/2 



9 = 9N- 



(1 + ^l+4ag/|gjvP) 



V2 



Equation |1] is much more tractable since gi^ may be 
calculated from the Newtonian potential as usual and then 
simply modified to give the correct MONDian acceleration 
at a given point in the field. 

Our galaxy is clearly neither planar, spherical nor cylin- 
drical and so the applicability of equation 21 may rightly 



be questioned. However, iBrada fc MilgrornI (^9^ demon- 
strated that equation may be used exactly for infinitesi- 
mally thin Kuzm in discs with Newtonian potential given by 
iBinnev fc Tremaine. 1987,1 : 



-GM 



^R^ + {a+\z\Y 



(5) 



Where G is the gravitational constant, a is the disc scale 
length and M is the mass of the disc. The reason that the 
Kuzmin potential can be used exactly is because it is an 
extremely special potential for which |V$jv| = /("I>jv)^. No- 
tice that in MOND, the force field is still the gradient of a 
scalar potential and so the MONDian field must be conser- 
vative; that is: V X p = 0. Thus a MONDian field may be 
generated via equation |21 from a Newtonian field provided 
that the Newtonian field satisfies the following constraint: 
V|V"I>jv| X V"I?]v — 0. This is satisfied exactly by the Kuzmin 
disc. 

In this paper we use the Kuzmin potential to study 
orbits in axisymmetric potentials in MOND. We compare 
these orbits to similar orbits in standard Newtonian me- 
chanics (the GDM model) using the Kuzmin potential plus a 
flattened spherica l logarithmic potential (to model the dark 
matter) given by llBinnev fc TremainellT987i) : 



$l{R, z) = —Vq In i + -{ — 7 I + constant 



(6) 



Where Rc is the scale length, 0.7 ^5 g ^5 1 is the halo flat- 
tening and 170 is the asymptotic value of the circular speed 
of test particles at large radii in the halo. 

We will also com pare these to th e more reali stic Galac- 
tic po tential used bv | Johnston et all (120051) and lLaw et alJ 



(|2QQ^. They use a logarithmic ha lo (equation Et, a 
Miya moto-Nagai potential for the disc iBinnev fc Trmnaine 



|l987l') and Hernquist potential for the bulge~ (lHernmtls1 
[1990): 



-GM 



bulge 



y/R^ + {a + Vz^ + b2)2 
r -I- c 



(7) 



(8) 



where a is the disc scale length as in equation El 6 is the 
disc scale height, M is the disc mass, c is the bulge scale 
length and Mbuige is the bulge mass. Notice that for b —> 
equation |71 reduces to the Kuzmin disc in equation El 



3 INITIAL CONDITIONS AND ORBIT 
SOLVING 

The mass distribution we use in MOND is the flattened 
Kuzmin disc (see equation EJ- For the CDM model we use 
a Kuzmin disc plus a logarithmic halo (see equation EJ- 
We present three CDM models: one with no halo flattening 
(CDM), one with q = 0.95 (f095CDM) and one with q = 0.9 
(f09CDM). We also compare these with th e bes t fit Milky 
Way p otential from Ijohnston et alJ (|2005^ and iLaw et all 
l|2005h (L05). The parameters used in all five models are 



For the Kuzmin potential, |V#]v| 
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Figure 1. Rotation curves for the MOND model (blue line), 
the CDM models (red line) and the LOS model (green line). 
The rotation curve which the Kuzmin disc would give with- 
out MOND is also overplotted for comparison (magenta line). 
The black data points show the mean of HI measurements 
of the rotation curve taken from [B ^ni^^ Lockma n 1 19^ ). 
IWeaver fc Williams7l974'). ' Weaver fc Williaml [l973ft . lMalhotral 
^1995^ and iKerr et al.i >198a) . 



given in table Inland are chosen to match the measured ro- 
tation curve of the Milky Way. 

Figure0shows the rotation curve for the MOND model 
(blue line), the CDM models^ (red line) and the L05 model 
(green line). The rotation curve which the Kuzmin disc 
would give without MOND is also overplotted for com- 
parison (magenta line). The black data points show the 
mea n of HI measurements of the rotation curve taken 
fromlBania fc LockrnanI lll^ 4). We aver fc WiUia-mj Jl974). 
IWeaver fc Wilhamj jlOTaTlMalhotral ^1995^ and lKerr et alJ 
It is important to note that we are not attempting 
to form an accurate model of the Milky Way in this paper, 
but instead we wish only to compare orbits in MOND and 
CDM galaxies. The Kuzmin potential is not the most ac- 
curate model for t he stellar d istrib ution of the Milky Way 
(see e.g. lCaldweiT fc Ostrike] lllflSJ^ ) and in both the CDM 
and MOND models, we have made no attempt to model 
the stellar bulge and bar although it is well known that 
they contribute signif i cantly to t he potential of the galaxy 
JCaldwell fc Ostrikeil Jl98li) and iDwek eraJl (119951) ). This 
can be seen in the difference between the L05 rotation curve 
with a bulge component (green line) and all of the other 
models. Interior to ~ lOkpc, all of the rotation curves devi- 
ate quite strongly from L05. However, since most Milky Way 
satellites orbit well outside of lOkpc, a potential model for 



the Milky Way which is accurate beyond this point should 
suffice. 

The equation of motion in MOND and in CDM is given 

by: 



(9) 



Where in MOND V"!? = g is calculated from equation |1| 

Equation |3| represents a set of coupled differential equa- 
tions which we solve n umerically using the fourth-order 
Runge-Kutta technique jPress et alJll993) with a timestep 
of 0.15Myrs. Reducing the time step was found to produce 
converged results, while for purely spherical potentials, the 
code was found to conserve energy and angular momentum 
to machine accuracy (better than 1 part in 10^)^. 



4 RESULTS 

We modelled the Sagittari us dwarf o r bit by fitting to the 
best-fit orbit presented in iLaw et all (|2005|) and ensuring 
that the current position and velocity of the dwarf matched 
observational constraints. This gives a current phase space 
position and velocity of the dwarf (in Galacto-centric right- 
handed coordinates) of x = 16.2, y = 2.3, z = —5.9,vx = 
238, vy — —42, vz — 222, where numbers are quoted in units 
of kpc and km/s respectively. This satellite phase space po- 
sition was then integrated backwards IGyr in time to match 
the trailing arm of the Sagittarius dwarf and forwards IGy r 
in time to match the leading arm. Unhke lLaw et al.1 J2005I1 . 
we did not allow the final phase space position of the dwarf 
to vary within observational constraints, but held this fixed. 
This is because we wish to measure the difference in orbital 
precession between the models, which is easier to do if the 
initial phase space coordinates are identical. 

Figure |5| shows orbital projections for all five models 
as shown in the legend, integrated over 2Gyrs. The galactic 
plane (not marked) is perpendicular to the z-axis. The top 
panels are for the best-fit Sagittarius dwarf orbit which is on 
a near-polar orbit around the galaxy jLaw et al ]l2005l) . The 
bottom panels are for a small-pericentre orbit which samples 
a wide range of the Galactic potential. The position of the 
Sagittarius dwarf now is marked on the top panels with a 
red square. The orbital pole precession is shown in the right- 
most panels. This measures the difference in angle between 
the vector perpendicular to the satellite's orbit initially and 
at a given time. 

All of the models produced very similar orbits for the 
Sagittarius dwarf and the differences are best seen in the 
right-most panel plots of the orbital pole precessions. The 
pole precession is a useful quantity to measure for the orbits 
because it is a strong function of how fiattened a potential 
is - this is why the differences between the orbits in each 
of the models shows up so strongly in the plot of the pole 
precession, whereas it is much harder to detect in the plots 
of the orbital projections. Small differences in the orbits be- 
tween models may be accounted for by altering the details 
of the visible component of the Milky Way potential (recall 



^ All of the CDM models will produce the same rotation curve 
since the force from the Galaxy on the satellite in the plane of 
the Galaxy is independent of the dark matter halo flattening, q. 



^ The orbits in the axisymmetric potentials used in this paper 
also conserved energy and the z-component of the angular mo- 
mentum to m achine accuracy - as expe cted for potentials with 
axisymmetry iBinnev fc Tremainelll987^ . 
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Table 1. Initial conditions. 
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Figure 2. Orbital projections for five models as shown in the legend, integrated over 2Gyrs. The galactic plane (not marked) is 

S ierpendicular to the z-axis. The top panels are for the best-fit Sagittarius dwarf orbit which is on a near-polar orbit around the galaxy 
Law et all2005! ). The bottom panels are for a small-pericentre orbit which samples a wide range of the Galactic potential. The position 
of the Sagittarius dwarf now is marked on the top panels with a red square. The orbital pole precession is shown in the right-most panels. 
For the Sagittarius dwarf orbit, the time marked is relative to its current phase space position. The observed trailing arms then trace 
out the orbital path the dwarf took over the past Gyr (hence the negative time), while the leading arms show what its path will be over 
the next Gyr (hence the positive time). 



that the Kuzmin disc used for the MOND model is only 
an approximation to the true potential of the Milky Way 
disc and bulge) ; or by altering the final phase space position 
of the Sagit tarius dwar f within observational constraints as 
was done bv lLaw et aP (2005) . The pole precession, however, 
may only be reproduced by changing how flat the potential 
is. In the CDM models, this may be achieved quite easily 
by using a more oblate dark matter halo. In MOND there 
is less freedom to do this. While changing the mass of the 
Milky Way disc can produce more or less precession, there 
are strong limits from stellar population models and from 
the rotation curve of the Milky Way as to how much this 
can be done. If MOND produces far too much, or far too 
little precession as compared with the Sagittarius stream, 
then we can rule it out as a viable alternative theory to 
dark matter. 

The CDM model (magenta line) produced the least pre- 



cession as can be expected for a near-spherical potential 
(recall that this model still contains a massive disc and so 
we should still expect some precession). The f09CDM and 
LOS models produced very similar results. This is because, 
with an orbital pericentre greater than lOkpc, the Sagittar- 
ius dwarf is not sampling a region of the Milky Way potential 
where the presence of the bulge is significant. The key point 
is that, surprisingly, the MOND model (blue line) produced 
near-idenUcal pole precession to both the LOS and f09CDM 
models - consistent with the best fitting orbit for the Sagit- 
tarius dwarf debris. If anything, the MOND model produced 
too little precession for the leading arm of the Sagittarius 
dwarf debris. This result seems surprising since in MOND 
all of the gravity is coming from the disc. We will discuss 
this further in section |S] 

There could, then, be an early indication that MOND is 
inconsistent with the Sagittarius stream orbit. However, cur- 
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rent data from the Sagittarius dwarf is only good enough to 
constrain the pole precession o ver the period -0 . 4 to .6Gyrs 
with an error of 2-3 degrees (|,Tohnston et al l l2nn,^ . This 
places a flattening oi q = 0.9 — 0.95 within the l.Scr error 
bars. Yet the difference in pole precession between q = 0.9 
and q — 0.95 is much larger than the difference between 
MOND and the f09CDM or L05 models. Oyer th e range - 
0.4 to 0.6Gyrs discussed in Ijohnston et al.l (|20ol), MOND 
produces near-identical results to the f09CDM and L05 mod- 
els. 

There is still some lively debate about what the best 
fitting orbit for the Sagittarius dwarf actually is (see section 
14. fl . However, it seems unlikely that it would be possible 
to produce more precession in MOND than that from an 
infinitesimally thin disc (the potential used in this paper). 
If this is true, then it may be possible in the future to rule 
out MOND on the grounds that it does not produce enough 
precession to match the Sagittarius stream data. 

The bottom panels show a more extreme orbit. Notice 
that now the orbit for the L05 model strongly deviates from 
the others. This is because this satellite orbit has a peri- 
centre of ~ 2kpc and so the satellite is now sampling the 
region of the Milky Way potential where the bulge makes 
a difference. The rest of the models, which have no bulge 
component, once again show very similar orbits. The pole 
precession in MOND is again best matched by the f09CDM 
model, while the f095CDM and CDM models show less pre- 
cession, as is expected from their rounder halo potentials. 

The MOND potential produces orbits which are very 
similar to those in a CDM halo with flattening of g ~ 0.9. 
This seems to be the case even for orbits with unrealistically 
small pericentres. This suggests that, for the Milky Way at 
least, it will be difficult to distinguish between MOND and 
a dark matter model for the Galaxy using satellite streams. 
Even if better data could be obtained for the globular cluster 
orbits (which orbit closer to the Galaxy), MOND can be 
expected to produce similar results to a mildly oblate dark 
matter halo. 

4.1 The Sagittarius leading arm velocity data 

There has been some debate in the literature over the veloc- 
ity data for the Sagittarius dwarf leading arm. We have so 
far in this paper referred only to the spatial data from the 
Sag ittarius stream and not the velocity data. 

iHelmJ (|20o3) has recently pointed out that the lead- 
ing arm velocity data are inconsisten t with a halo flattening 
of q = 0.9 advocated by Johnston et al.l l|200^ and sug- 
gest that a prolate halo (g ~ 1.25) would provide a bet- 
ter fit if the velocity d ata are taken into account. However, 
Ijohnston et al.l (12005^ argue that a prolate halo leads to an 
incorrect value for the orbital pole precession; in fact a pro- 
late halo causes the orbit to precess in the opposite direction 
to that observed. They point out that it is difficult to obtain 
the correct amount of pole precession without altering the 
underlying potential (as can be seen in the pole precession 
plots in figure ; whereas one could conceivably alter the 
velocities of the stars in the plane of the strea m through 
secon d order effects such as dynamical friction jLaw et alJ 
[200^. Could MOND perhaps reconcile the discrepant lead- 
ing arm velocity data? 

Figurel^shows the the line of sight velocity of the Sagit- 



tarius stream as viewed from the sun as a function of its 
longitude on the sky. The colours for the models are as pre- 
viously: MOND is blue, LOS is green and f09CDM is red. 
The extra line shown is for a prolate CDM model (yellow 
line) with q = 1.25. It is the yellow line which provides the 
best fit to the velocity data for the Sgr stream, particularly 
the leading arm data (left most curves on the plot). Notice 
that the MOND model (blue line) agrees well with the L05 
and f09CDM models as before (red and green lines), but 
that all of the models differ greatly from the prolate model. 

MOND does not solve the problem of the discrepant 
leading arm velocity data for the Sagittarius stream. As with 
the spatial data for the stream, MOND produces a near- 
identical orbit to the L05 and f09CDM models. If it can be 
shown that the Milky Way halo (or a ny other galax y halo) 
must be prolate, this, as pointed out bv lHelml (|200^, would 
be difficult to reconcile with MOND. 



5 DISCUSSION 

5.1 Model assumptions 

Perhaps the biggest assumption in this work is the highly 
specific choice of potential for the galaxy disc. While this 
could be problematic for a detailed study of satellite orbits 
in the Milky Way, we have implicitly considered the case of 
maximal precession in MOND. It would be difficult to imag- 
ine a more axisymmetric potential than an infinitesimally 
thin disc. This does mean, then, that should a dark matter 
halo be discovered which is significantly more oblate than 
q — 0.9, this would be difficult to reconcile with MOND. 

We have also neglected dynamical friction and not per- 
formed a detailed N-body simulation of the stripping of stars 
from the Sagittarius dwarf galaxy. As such, this work should 
not be taken as conclusive evidence that the Sagittarius 
stream is consistent with MOND. 

Finally, we should note that MOND is probably the 
bes t-studied but not th e only alternative gravity theory (see 
e.g. iDrummondl (|20o3) and [Moffat ( 2004)). Some of these 
theories predict a Keplerian fall off at large radii in the rota- 
tion curve, similar to CDM models (see e.g. Moffat (2004)). 
These may be even more difficult to rule out using tidal 
streams. 



5.2 Why does such a flat mass distribution in 
MOND produce such a round potential? 

This point has been discussed in some detail in iMilgrorJ 
( 2001), but is perhaps best illustrated by direct integration 
of equation 0] for a Kuzmin disc. In the deep-MOND limit, 
l-i(x) X and we find from equations |1] and 13 

^,.^M^^ln{R'^{\z\+af) (10) 

Which is very nearly identical to the ffattened logarith- 
mic halo (see equation |SJ. Thus, highly ffattened mass dis- 
tributions in MOND do not produce highly ffattened po- 
tentials. In fact, it is the spherical nature of the Kuzmin 
potential in the deep-MOND limit which leads to MOND 
producing slightly too little precession in the Sagittarius or- 
bit (see figurel^J. From figure^ we can see that the MOND 
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Figure 3. Line of sight velocity of the Sagittarius stream as viewed from the sun as a function of its longitude on the sky. The colours for 
the models are as previously: MOND is blue, LOS is green and f09CDM is red. The extra line shown is for a prolate CDM model (yellow 
line) with q = 1.25. It is the yellow line which provides the best fit to the velocity data for the Sgr stream, particularly the leading arm 
data (left most curves on the plot). Notice that the MOND model (blue line) agrees well with the L05 and f09CDM models as before 
(red and green lines), but that all of the models differ greatly from the prolate model. 



rotation curve (blue line) is flat beyond ~ 20kpc, indicating 
that it is then in the deep-MOND limit. The other rota- 
tion curves for the CDM models are all falling at these radii 
rather than fiat - a property which cannot be achieved in 
MOND, since it is a theory set up to produce flat rather 
than falling rotation curves at large radii. Thus interior to 
~ 20kpc, all of the models agree quite well in their rota- 
tion curves and the corresponding orbit for the Sagittarius 
trailing arm looks very similar. This can be seen in the good 
agreement from -1 to OGyrs in the orbital pole precessions 
(figure 121 top right panel). However, the leading arm orbit 
- which can be seen in the pole precession plot from to 
IGyrs - does not agree so well. For this part of the orbit, the 
Sagittarius dwarf moves out towards apocentre and samples 
the region of the potential where the Milky Way is fully in 
the MOND regime and where the rotation curve (in MOND) 
is flat and near-spherical. 

5.3 Exploiting flattened eUiptical galaxy 
potentials 

iBuote fc Canizarej Jl994) and lBuote et al.l ll2002l) have re- 
cently shown that the observed flattening of hot X-ray gas in 
elliptical galaxies may also be used to place tight constraints 
on MOND. They argue that, if the hot gas in NGC720 is in 
hydrostatic equilibrium, then Vpgas = — pV$ which implies 
that Vp X V$ = 0. Thus, the X-ray isophotes from the hot 
gas in NGC720 trace the gas density, which in turn traces 
the underlying gravitational potential. By deprojecting the 
stellar potential they show that the stars in NGC720 cannot 



produce a flat enough potential in MOND to produce the 
observed X-ray isophotes. 

While they have to assume hydrostatic equilibrium and 
some simple form for the deprojected stellar potential, their 
assumptions are quite conservative. They flnd that MOND 
produces potentials which are too spherical at large radii - 
similar to what may be the case here for the leading arm 
Sagittarius dwarf debris. 



5.4 What are the prospects for constraining 
MOND using tidal streams? 

Cold dark matter halos, as modelled in N-body cosmologi- 
cal simulations, are typically 2:1 triaxial systems. This would 
make the Milky Wa y with q ~ 0.9 qui t e rare . However re- 
cent simulations by iKazantzidis et all ll2004l) showed that 
the dissipation of baryons changes the inner galactic ha- 
los to be nearly spherical and consistent with the flattening 
predicted from the Sagittarius stream. The amount of flat- 
tening depends sensitively on the fraction of baryons that 
undergoes slow dissipation to form the galactic disk. 

MOND mimics halos with q ~ 0.9, while CDM produces 
similar halos as a result of gas cooling and galaxy formation. 
This will make it difficult to differentiate between MOND 
and CDM theories using halo flattening, even with a large 
statistical sample of halo shapes. 
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6 CONCLUSIONS 

Wc have compared orbits in a thin axisymmetric disc poten- 
tial in MOND to those in a thin disc plus near-spherical dark 
matter halo predicted by ACDM cosmology. We demon- 
strated that the amount of orbital precession in MOND 
is very nearly identical to a similar CDM galaxy with a 
logarithmic-halo with flattening q=0.9, consistent with re- 
cent constraints from the Sagittarius stream. Since very flat- 
tened mass distributions in MOND produce more spheroidal 
potentials than in standard Newtonian mechanics, we have 
shown that it will be very difficult to use the tidal debris 
from streams to distinguish between a MOND galax;y and a 
standard CDM galaxy with a mildly oblate halo. 

If a galaxy can be found with either a prolate halo, or 
one which is more oblate than q ~ 0.9 this would rule out 
MOND as a viable theory. Improved data from the leading 
arm of the Sagittarius dwarf - which samples the Galactic 
potential at large radii - could rule out MOND if the orbital 
pole precession can be determined to an accuracy of the 
order of ±1°. 
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